Reciprocity of Networks with Degree Correlations and Arbitrary Degree Sequences 
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Although most of the real networks contain a mixture of directed and bidirectional (reciprocal) 
connections, the reciprocity r has received little attention as a subject of theoretical understanding. 
We study the expected reciprocity of networks with an arbitrary input and output degree sequences 
and given 2-node degree correlations by means of statistical ensemble approach. We demonstrate 
that degree correlations are crucial to understand the reciprocity in real networks and a hierarchy of 
correlation contributions to r is revealed. Numerical experiments using novel network randomization 
methods show very good agreement to our analytical estimations. 
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I. INTRODUCTION 

Most of real networks combine both unidirectional and 
bidirectional (reciprocal) connections. This directed na- 
ture is often obviated, e.g. networks are symmetrized 
for algorithmic convenience and network models largely 
ignore the directionality for analytical simplicity. How- 
ever, the directionality is known to be relevant, e.g., ro- 
bustness against environmental changes of metabolic net- 
works seems to arise from evolutionary pressure on the di- 
rections and weights of the metabolic fluxes [l[ . The for- 
mation of functional communities and hierarchies in the 
cerebral cortex is mediated by the presence of reciprocal 
and unidirectional connections 0. The dynamical sta- 
bility in complex networks, e.g. ecological systems 0, Q 
and synchronization of coupled oscillators, is commonly 
assessed by the eigenvalue space of the Jacobian or Lapla- 
cian matrices [Ij. When networks are directed these ma- 
trices will have complex eigenvalues, what influences both 
the stability and the dynamical organization far from the 
equilibrium state. 

The network reciprocity r is classically defined as 
r = ^j— where is the number of directed links 
.s t that also have a reciprocal (bidirectional) coun- 
terpart s <— i, and L is the total number of directed 
links. In networks without self-loops, reciprocal links 
form the cycles of lowest order and are therefore, im- 
portant as a natural measure of feedback in the network. 
Recently, r of the Wikipedia networks [2^1 (for differ- 
ent languages) was found to be very stable over a wide 
range of network sizes Q which signals its relevance for 
the structure or functionality of the networks. In Q it 
was shown that reciprocal connections carry most of the 
topological information of the WWW. The formation of 
the giant-component in directed networks is facilitated by 
reciprocal connections Q. Despite the extensive model- 
ing efforts during the recent years to reproduce realistic 
features of networks, models have largely ignored reci- 
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procity. Only in [To| a general class of random networks 
with prescribed r has been presented. 

When analyzing real networks it is important to test 
whether measured values are significant or not. Typi- 
cally, real network measures are compared to the prop- 
erties of complete random networks of the same size N 
and number of links L. However the degree distribution 
of most real networks largely differs from the aforemen- 
tioned random networks and more reasonable compari- 
son is desired. Derivation of analytical expressions for 
expected measures under conditions of arbitrary degree 
sequence is difficult. Further assumptions are usually 
introduced, e.g. scale-free [lH or exponential degree dis- 
tributions [I4I . As real networks do never exactly belong 
to a model class, if even close, for real applications a nu- 
merical approach is the only solution in most situations. 
Ensembles of maximally random networks can be gener- 
ated with the same input and output degree sequences 
and the ensemble average properties can be calculated. 
Unfortunately, generation of such ensembles is compu- 
tationally very demanding for large networks. In this 
paper, we present analytical expressions for the expected 
reciprocity of directed networks that can be evaluated 
using only information measurable from the specific real 
network under study, and thus, overcome the problems 
discussed above. 

In the context of social networks, several redefinitions 
of r have been introduced that account for biases of 
the experimental conditions 0, [l^. Recently, a simi- 
lar redefinition based on correlations of adjacency ma- 
trix and its transpose, was presented It evaluates 
the reciprocity with respect to the density of connections 
a = L/{N{N - 1)). This results from the fact that for 
a maximally random network with N nodes and L links, 
r equals a. However, the specific degree sequence is ex- 
pected to affect the number of reciprocal links, i.e., a node 
with both large input degree ki and large output degree 
ko has a higher tendency to form reciprocal links. This 
is expressed by the 1-node degree correlations between ki 
and ko of individual nodes. Similarly, as r involves the 
pairwise connectivity of nodes, the correlation between 
the degrees of neighboring nodes, 2-node degree corre- 
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FIG. 1: 2-node 2-degree correlations {2n2d) of neighboring 
nodes in directed networks. Links corresponding to correlated 
degrees are colored black. 



lations should also be relevant. Imagine a link s — > t 
connecting a source node s with input and output de- 
grees {ki,ko) to a target node t with degrees {qi, Qo). The 
2-node degree correlations exist when any of the {ki, ko) 
degrees of node s is correlated with any of the {qi, qo) de- 
grees of node t. When all four values arc correlated, then 
both 1-node {In) correlations and the 2-node 4-degree 
{2nAd) correlations are present. The class of 2-node 2- 
degree {2n2d) correlations is depicted in Fig. [T] 

In this paper we study the expected reciprocity (r) of 
networks with prescribed degree sequences and arbitrary 
2-node degree correlations. We consider complex net- 
works as members of the statistical ensemble with given 
node degree sequence and degree correlations, and we 
calculate the expected reciprocity of such ensembles in 
the thermodynamical limit. We find that degree correla- 
tions explain almost completely the observed r of some 
real networks. In other examples, larger discrepancies 
indicate the presence of additional internal structure. 



II. GENERAL RESULTS 

In order to analytically estimate (r) under different 
correlation structures, we characterize the real directed 
networks by: the number of nodes N, the number of 
links L, the number of nodes N{ki,ko) = N{k) having 
in-degree ki and out-degree ko, and the number of di- 
rected links L{k — *■ q) pointing from nodes with degrees 
{ki,ko) to nodes with degrees {qi,qo). All these proper- 
ties are easy to measure in a real directed network and 
contain all the relevant information about the degree cor- 
relations. We use frequencies of these properties as their 
probabilities and calculate the expected number of re- 
ciprocal links {L^). Remind that by definition, (r) is 
related to {L^) by (r) = (L^) /L. 

Under the class of 1-node and 2-node degree correla- 
tions here assumed, a network is considered as maximally 
random when any of the nodes with degrees k is equally 
likely connected to any of the nodes with degrees q. If 
a network contains L{k q) such links, the probabil- 
ity that any of them connects randomly chosen nodes of 
degree k and q respectively is, 

, L{k ^ q) 



p{k 



N{k)N{q) 



(1) 



in the thermodynamical limit |27| . The denominator 
N{k)N{q) is the number of all possible connections be- 
tween nodes with degrees k and nodes with degrees q. 

Again, if the network has L{k — > q) links of the type 
k ~> q, the expected number of reciprocal k q links 
is then {L{k ^ q)) = L{k q)p{k ^ q). The over- 
all expected reciprocity ri„2n of the whole network is 
obtained by summing {L{k <-> q)) over all k, q degree 
combinations: 

I y-^L{k-^ q)L{k ^ q) 
J. 2^ 



k,q 



N{k)N{q) 



(2) 



Note that in general L{k q) L{k ^ q). Taking 
frequencies of nodes P{k) = N{k)/N and frequencies of 
links 'P{k q) = L{k q)/L as probabilities in the 
thermodynamical limit, Eq. ([2|) reads, 

L ^ V{k^q)r{k^q) ^ 



k,q 



P{k)P{q) 



The contribution of the correlation structure is accounted 
by the sum, V{k q) accounts for both 1-node and 2- 
nodc correlations, and P{k) only for the 1-node correla- 
tions. When all four degrees are independent V{k 
q) = P{h)koP{ko)q^P{qr)P{qo)l {kf where {k) = 
L/N is the average degree and P{k) = P{ki) P{ko). 
Then reciprocity reduces to the density of links L/N'^, 
i.e., the expected reciprocity of uncorrelated random net- 
works. 



III. SPECIAL CORRELATION CLASSES 

Equations ([2]) and ([3]) are general formulas that ac- 
count for all 1-nodc and 2-node degree correlations. 
These equations can be reduced to consider only desired 
special classes of correlations and thus explore the contri- 
bution of individual correlation types to r. In this section 
we present detailed derivations for all 8 possible combi- 
nations of 1-node and 2-node correlations and in Table U 
the main results are summarized. Note that in Table U 
only Tc, the contribution of the correlation structure is 
shown. The expected reciprocities arc obtained by mul- 



tiplying with the density of links: (r) 
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this section, the usual product rule of joint probabilities 
in terms of conditional probabilities will be used, 

P(Xi,X2,...,X„) = P(Xi) X P{X2\X{) X ... 

... X P{Xn\Xi, X2, Xn-l). 

With this rule in mind, the joint probabilities associated 
with the link statistics 'P{k —> q) = p(ki, ko, qi, qo) can 
be expressed, for example, as: 

r{k ^ q) = V{ko qi)P{h I ko -> qi)P{qo \ ko ^ q^, h) 



A. 2-node out-/in-degree correlations 

Suppose that a network has significant 1-node corre- 
lations and 2-node correlations only between the out- 
degree ko of source nodes and the in-degrce qi of the 
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target nodes, see Fig. [Tfa), while other possible 2-nodc 
correlations are negligible. In this case, the probabilities 
can be approximated by: 

Vik ^ q) « V{ko ^ Qi) P[h\ko) PiQoh) 
V{k ^ q) « r{K ^ qi) P{ko\h) P{qMo) 

Using conditional probabilities for the degrees of individ- 
ual nodes, e.g. P{ki\ko) = P{k)/P{ko), Eq. ^ reduces 
to: 

_J^Sr ^(^° ^ ^ qo)P{k)P{q) 



1-nodc 2-node 



k,q 



P{h)P{ko)P{q^)P{qo) 



(4) 

Additionally, if the 1-node correlations are negligible, the 
degrees of individual nodes become independent: P{k) — 
P{h,ko) P{ki)P{ko). Then Eq. gj becomes: 

r2n:o/i = — ^(fco ^ qi)V {h ^ qo) = —, (5) 



N2 

k,q 



that equals the density of connections a in the thcrmo- 
dynaniical limit. This means that the 2-node out-/in- 
degree correlations do not contribute to reciprocity in 
the absence of 1-node correlations. This is not a gen- 
eral case since other classes of 2-nodc correlations largely 
contribute to r. 



B. 1-node degree correlations 

Starting from Eq. ^ , we can alternatively remove the 
remaining 2-node in-/out-correlations and obtain a gen- 
eral expression for the expected reciprocity ri„ due to 
the 1-node correlations alone. In this case, the num- 
ber of all possible output connections from source nodes 
with degree ko is ko N{ko) and the number of all possi- 
ble input connections to target nodes with in-degrec qi is 
qiN{qi). Hence, the probability that one link connects 
a node with out-degree ko to a node with in-degree qi 
is V{ko q.) - SMS^ ^ k.Pik^),.PM ^ 

tion ([4]) reduces to: 

L (kiko) 



rir. 



Ar2 



{ky 



(6) 



Here r, a two node property, is determined by the single 
node characteristics arising from the specific in- and out- 
degree sequences. We remind that, in the literature, it 
is common to randomize networks by methods that con- 
serve the degree sequences in order to obtain expected 
values accounting for the real degree distribution. Equa- 
tion ^ is of relevance for significance testing because it 
is the theoretical estimation of the expected reciprocity 
in such a typical case. 



C. 2-node out-/out-degree correlations 

Following a similar approximation we can calculate the 
expected reciprocity r2n:o/o due to the 2-node correla- 
tions, Fig. [Ijb), between the output degree ko of the 
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k,,q, P{k,)P{qi) 
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yes 
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k,q 

V{ki^qo)V{ka<-qi)P(ki,ka)P{qi.qa) b 
kj^kp ^oM^ ^i,Qo ^0:9i 



"The averaging in the formulae is performed over in degrees of 
source nodes and the out-degree of the target nodes. 

'^o fc = Sfc' k'oPi^'o^ ki) and similar for all other averages of this 
type. 

TABLE I: Reduced formulae for the expected reciprocity of 
different combinations of 1-node and 2-node correlations. The 
expected reciprocity is (r) = Tc- Note that the result for 
the 2-node out-out and in-in correlations are independent of 
the 1-node correlations. 



source node and the output degree qo of the target node. 
However, in this case a few steps need to be carefully 
considered. We rewrite the link probabilities as: 

P(k ^ q) = V{ko ^ qo)P{ki\ko — > qo)P{qi\ki, 

The term P{qi\ki, ko — > qo) introduces dependence of in- 
dcgrce of the target node qi on the in-degree of the source 
node ki. We are assuming that such correlations are 
negligible and therefore this term can be rewritten as: 
P{qi\ko qo)- The second term P{ki\ko — > qo) can be 
written as P{ki\ko), because the in-degree of the target 
node qi and the in-degree of the source node ki are, again, 
not correlated. To proceed, it is necessary to use the fact 
that we are calculating expectations on graphs, and that 
the final expression should has a form which can be cal- 
culated using only the assumed knowledge of the network 
structure, i.e. the frequencies of nodes P{k) = N{k)/N 
and the frequencies of links 'P{k q) = L{k — > q)/L. 
Following this line of reasoning the third , already ap- 
proximated term, should be carefully rewritten. What is 
the probability that a node will have in-degree qi given 
that (i) it has an input link and (ii) it has an out-degree 
Qo? Without the information (i) this probability is sim- 
ply P{qi\qo), but we have information that such a link 
does exist. The probability that, following a link between 
nodes with degrees ko and qo, it will run into a node with 
in-degree qi is just the number of links i(— *■ q) which en- 
ter all nodes with degrees q divided by the number of all 
links i(— > qo) that enter the nodes with out degree qo- 
These numbers of links can be expressed as qiN{qi,qo) 
and Yl,q' q'i^iq'i^qo) respectively. Thus, 

qiN{qi,qo) qiP{qi,qo) 



P{qi\ko 



qo) 



Y.,'q[N{q[.qo) E^'q^m^qo) 



Substituting all these approximations in Eq. ([3]), 
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Viko ^ qo)V{K ^ qo)hP{k)q,P{(i) 



r2n:o/o ^ k[P{k[, K)q[P{q[, qo)P{K)P{qo) 

Finally, the summation terms over in-degrees in the nu- 
merator and the denominator cancel out and we obtain a 
final expression for the expected reciprocity under 2- node 
out-/out-dcgree correlations: 



L 



E 



V{ko^qo)V{ko^qo) 
P{ko)P{qo) 



(7) 



Interestingly, this expression is independent of the 1- 
node correlations even if we did not explicitly assume 
it. The expression for 2-node in-/in-degree correlations 
conserved, Fig. [T]Jc), is the same only with ko and qo 
replaced by ki and qi, see Table HI 



D. 2-node in-/out-degree correlations 

We now describe the case in which significant in-out 
correlations, Fig.[T]^d); arc present in the network. Such 
correlations are supposed to influence reciprocity consid- 
erably because the probability that a given link of the 
type k ~^ q has a reciprocal counterpart fc ^ q is di- 
rectly proportional to the in-degree ki of the source node 
and the out-degree qo of the target node. Following the 
previous line of reasoning, the link probability can be 
expressed in this case as: 



V{\<i ^ q) = V{h ^ qo)P{ko I k, ^ qo)P{q^ \ h ^ 

koP{ko-! ki) 



and the conditional terms as: 



P{ko I ki 



P{q^ I h 



Y.k'^k',N{kM T.k-yoP{k'o,ki) 
koP{ko^ ki) 



ko,ki 

qiN{qo,qi) 



qiP{qo,qi) 



T.,'AN{qo,q[) T.,'Anqo.q[) 

qiP{qo,qz) 



Note that the expression ko^kt = J2k' k'oP{ko, h) is not 
the average out-degree of nodes with in-degree ki, be- 
cause P{k'^, ki) is the joint probability and not the condi- 
tional probability. The average out-degree of nodes with 



in-degree ki is {ko)f.. 



P(k.) 



Using the above rela- 



tionships we obtain from Eq. ([3]) an expression for the 
expected reciprocity ri^^n-.i/o = "^i/o due to 1-node and 
2-node in-/out-degree correlations: 

L ^ hkoq^qoV{h ^ qo)V{ko ^ q,)^(k)-P(q) 



k,q 



i,ko ^o,ki Qi^Qo Q.o,qi 



.(8) 



Additionally, if the 1-node correlations are removed, the 
expected reciprocity rm-Ajo due to the 2-nodc in-/out- 
degree correlations alone is: 



_ J^Sr^ kikoq,qoV{ki qo)'P{ko ^ q^) 
r2n./o - 2^ i^ko),^ ih),^ (qo)^^ 

Note that without 1-node correlations (fco)^ = 

j:,, k'^p{k'^,h)/p{ki)^{ko). 



IV. APPLICATION TO REAL NETWORKS 

The class of degree correlations chosen in this paper 
is not only very interesting from the theoretical point 
of view, in this section we pay attention to its practical 
relevance. First, our results are applied to several real 
networks and the impact of degree correlations on reci- 
procity is discussed. Finally, in order to proof the valid- 
ity of our equations as expectation values, the theoretical 
estimations are compared to the empirical ensemble av- 
erages of random networks. Therefore, we present several 
algorithms to generate random networks conditional on 
different correlation classes. 



A. The reciprocity of real networks 

All the analytical expressions summarized in Table |l] 
can be directly estimated by measuring the necessary 
statistics out of a real network. For different cases spe- 
cific quantities need to be counted, N{ki), N{ko), N{k) 
or L{ko qi), L{ki qo), L{k q), etc. Then, 
the frequencies like 'P{k ^ q) ~ L{k — > q)/L or 
P{ki) = N{ki)/N are introduced in the formulas of Ta- 
ble U to obtain the expected reciprocity under desired 
combination of 1-node and 2-node degree correlations. 

For each of the real networks in Table HIl we have cal- 
culated the expected reciprocities under different corre- 
lation classes, Eqs. ([3]), (O, (O and ([7]). For comparison, 
we also show the density of connections d, i.e., the reci- 
procity of complete random networks of the same size 
and number of links L. As observed in Table [TTl 
our general estimation of reciprocity ri„2n alone can al- 
most completely explain r of many networks, e.g. World- 
Trade- Webs and most of the food webs analyzed. It also 
makes a very good approximation for cortical and neural 
networks. However, the large discrepancy in the case of 
Wikipedia web-sites suggests the presence of additional 
internal structure in the network rather than 1-node and 
2-nodc degree correlations. For example, the C. elegans 
neural network, the cortical and the Wikipedia networks 
are known to have modular and hierarchical structure. 
In all these cases the real r is larger than the expected 
'i'in2n- Nevertheless, our results demonstrate that degree 
correlations are crucial to understand the reciprocity in 
real networks. In its closest approximation. World Trade 
Webs, ri„2n is roughly twice as large as d, the expected 
reciprocity of an equivalent complete random network. 
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Network r 








a 


World Trade Webs [15] 










Year 1948 0.823 


0.812 


0.768 


0.707 


0.382 


Year 2000 0.980 


0.958 


0.883 


0.813 


0.560 


Neural Networks 










C. Elegans [16] 0.433 


0.329 


0.071 


0.060 


0.033 


Cortical Networks 










Cat [IT] 0.734 


0.659 


0.473 


0.390 


0.300 


Macaque [IS] 0.750 


0.645 


0.287 


0.230 


0.155 


Food Webs [19| 










Little Rock lake 0.0339 


0.0323 


0.0365 


0.0501 


0.0743 


Grassland 0.0 


0.0 


0.0077 


0.0079 


0.0179 


St. Marks sea. 0.0 


0.0075 


0.0500 


0.0703 


0.0948 


St. Martin Isl. 0.0 


0.0016 


0.0419 


0.06765 


0.1131 


Silwood Park 0.0 


0.0 


0.0002 


0.0160 


0.0155 


Ythan estuary 0.0034 


0.0050 


0.0335 


0.0531 


0.0330 


Wikipedia Website [7] 










Spanish 0.3517 


0.1466 


0.0322 


0.0056 


0.0004 


Portuguese 0.3563 


0.1207 


0.0168 


0.0084 


0.0004 


Chinese 0.3668 


0.1556 


0.0256 


0.0096 


0.0010 



TABLE IL Measured reciprocity r of several real networks, 
and theoretically expected reciprocities due to different corre- 
lation structures, (i) 1-node and all 2-node correlations ri„2ii. 
(ii) 2-node out-out r2n:o/o (no 1-node). (iii) 1-node correla- 
tions Tin (no 2-node). And (iv) the density of connections 
a. 



In other cases ri„2n and a differ up to 3 orders of mag- 
nitude. 

The values of ri„2n are always followed by the con- 
tribution of 2-node out-out correlations r2n:o/o: Eq. ([7]). 
With the exception of Silwood Park and Ythan Estuary 
food webs, values of ri„, Eq. are all closer to the 
real r than the density of links a. These results reveal a 
hierarchy of the types of degree correlation with respect 
to their approximation level for explaining reciprocity. 



B. Numerical Corroboration 

Finally, we will proof that our analytical expressions 
are valid expectation values of the network reciprocity. 
For the real networks in Table HIl we have generated en- 
sembles of maximally random networks under different 
conditions of 1-node and 2-node degree correlations. The 
ensemble average reciprocities have been compared to the 
analytical results showing excellent agreement between 
experimental and theoretical expected values. Fig. [21 The 
generation of random networks with desired degrees and 
correlations is not trivial, but we have developed three 
novel rewiring / generation algorithms for that purpose. 
In the absence of analytical results for other graph mea- 
sures, e.g. clustering coefficient, average pathlength, etc., 
the following algorithms arc also useful to empirically cal- 



culate the expected value of any network measure under 
desired conditions of degree correlations. 



1. Random networks with desired degree sequences and 
2-node degree correlations, ri„2»i 

Given a real network of size N and L links, we can 
measure its degree distribution N{k) and the 2-node cor- 
relation structure L{k — > q). With this information in 
hand, it is possible to generate a maximally random net- 
work with such properties. To an initially empty net- 
work of size N , its nodes are randomly assigned their 
final degrees k following the distribution N{k). Then, 
links are introduced at random but following carefully 
considered steps. First, one source node s is chosen at 
random. We know that s has been assigned to have fi- 
nal degrees k' . Because of the 2-node correlations, s can 
only connect to nodes with particular degrees q' such 
that L{k' ^ q') > 0. A list of possible target nodes is 
constructed by taking only those nodes assigned to have 
final degrees q' where L{k' — > q') > 0. From this list one 
target node t is chosen at random and the link s — » <: is 
include to the initially empty network. Note that a node 
with input and output degrees k = (fc,;, ko) can only be 
eligible kg times as source and ki time as target, other- 
wise the distribution N{k) will not be conserved. Thus, 
if the quantities L{k q), and ki and ko of each node 
are adequately updated during the process, only L iter- 
ations are required to construct the maximally random 
network. Described as it is, this method allows for the 
introduction of self-loops and multiple- links. Avoiding 
them is far from trivial. 

To proof the validity of our general theoretical re- 
sult ri„2„, Eq. ([3]), ensembles of 100 random networks 
have been generated for each of the real networks in Ta- 
ble [IT] following the method described above. The results 
in Fig. [2Ia) show an excellent agreement between ri„2n 
(horizontal axis) and the empirical ensemble average re- 
ciprocities (ri„2Ti) (vertical axis). All the generated ran- 
dom realizations have been positively tested to have the 
same N{k) and L{k — » q) as the original real networks. 
From all the networks in Table HI] none but the C. Ele- 
gans contain self-loops and multiple- links. Therefore, the 
impact of the self-loops and multiple-links introduced by 
the random generation method was also tested. In the 
random realizations out of the Silwood Park food web, 
on average, only 4.5% of the L links formed self-loops 
or multiple- links. The realizations out of the Wikipedias 
contain on average less than 2% of such links, and the 
realizations out of the cortical networks and the world- 
trade- webs contain less that 1%. These small differences 
have very little impact as observed in the good agree- 
ment hown in Fig. [DJa) between our theoretical estima- 
tion ri„2n and the experimental measures (ri„2n)- 
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FIG. 2: Numerical corroboration of theoretical estimations. 
Experimentally measured reciprocities (vertical axes) and our 
theoretical estimations (horizontal axes), a) All 1-node and 
2-node correlations conserved, b) Only 2-node out-out cor- 
relations conserved, c) Only 2-node out-in correlations, d) 
1-node correlations. All data points are averages of 100 real- 
izations. Food webs (x), C. Elegans (•), WTW (T), cortical 
networks (a) and Wikipedias (*). Dashed lines are references 
of perfect coincidence. 



2. Rewiring method that conserves 2-node out-out 
correlations, r^n-.ojo 

In this case, we opted for a rewiring algorithm. Such 
methods start from a given real network and stepwise 
randomize its connections under specific rules that con- 
serve desired properties. In order to obtain maximally 
random networks that conserve uniquely the original 2- 
node out-out degree correlations, Fig. [ijb), we proceed 
in the following manner. One of the links in the network 
is chosen at random s t\. Nodes s and t\ have out- 
degrees fco and qo respectively. From all the nodes in the 
network with out-degree qo a new target node ti is ran- 
domly chosen. If all conditions to avoid the introduction 
of self-loops and multiple links are satisfied, the old link 
is destroyed s -^^ t\ and the new link is included s ^ t2- 

After several iterations, all 2-node degree correlations 
except for the out-out will be randomized because the 
rewiring step is blind to all other correlations. A rele- 
vant question arises when applying rewiring algorithms: 
how long should the process run so that resulting networks 
are maximally random? After some finite number of it- 
erations the network reaches a maximally random state 
and any successive rewiring will lead to an statistically 
equivalent random network. Once this state is reached, 
all network measures converge to an stable value. We 
have applied several levels of rewiring to the real net- 
works studied in this paper and the reciprocity has been 
measured at each level. As observed in Fig. [Slja) the 
average values of r reach a stable point. Any other net- 
work measure will follow the same behavior. Obviously, 
the number of necessary iterations is proportional to the 



number of links L and, for this particular method, all 
networks reach a maximally random state after 2L iter- 
ations. 

To proof the validity of our theoretical expression 
^2)i:o/oi Eq. ([7]), as expected reciprocity of networks with 
prescribed 2-node out-/out-degree correlations, we gen- 
erated ensembles of 100 rewired networks out of the real 
networks in Table [III For security, 4L rewiring itera- 
tion steps were used in all cases. The ensemble average 
reciprocities {r2n:o/o} were calculated. The comparison 
between our theoretical estimates and the empirical re- 
sults, Fig.[2)^b), shows again excellent agreement. 



3. Rewiring method that conserves 2-node out-in 
correlations, r2n:o/i 

In order to proof the validity of our theoretical r2n:o/ii 
we designed yet another rewiring algorithm that con- 
serves uniquely the original 2-node out-in correlations 
Fig. [Ija) of a real network. First, two nodes are se- 
lected at random, si and S2. These nodes are divided 
in two halves, one containing all in-links and the other 
containing all the out-links, i.e. si = s™ U 5°"* and 



S2 



U 



Then, the in-halves and the out-halves 



are switched forming two new nodes, S3 = U ^nd 
S4 = S2" U ^i"*- Cases that would introduce self- loops 
are carefully discarded. The resulting randomized net- 
works conserve the in-degree N{ki) and the out-degree 
N(ko) distributions, and the number of links L(ko — > qi) 
from the original network while the rest of correlations 
are randomized. 

The stability of the rewiring process has been tested in 
a similar manner as in the previous section and the re- 
sults are shown in Fig. [3I^b). In this case, each iteration 
step rewires several links so that the algorithm is much 
faster. After only 0.3L iterations all networks reach a 
maximally random state. The theoretical expected reci- 
procity r2n:o/ii Eq. ([5]), is compared to the empirical en- 
semble averages (r2n:o/i) of 100 rewired networks and 
the results shown in Fig.[2)^c). All networks were rewired 
QAL times. 



4-. Conserving degree sequences, ri„ 

Finally, a well known rewiring method was used to 
randomize all 2-node correlations while degree sequences 
N(k), 1-node correlations, are conserved [2(1. This 
method arises from earlier approaches [U, |22|| and con- 
sists in randomly choosing two links, Si — > ti and 
S2 t2, and exchanging them, si t2 and S2 — > ii. In 
Fig-HJd) the expected reciprocities (ri„) from the ensem- 
ble averages arc compared to our analytical estimations 
rm, Eq. dSI). 

All the results in this section proof that our analytical 
expressions are valid expected reciprocities under the cor- 
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FIG. 3: Stability of rewiring algorithms. With increasing 
number of randomizing iterations structural properties reach 
stable values, here r is shown, a) Algorithm that conserves 
the 2- node out-/out-degree correlations, b) Algorithm that 
conserves the 2-node out-/in-degree correlations. All data 
points are averages of 50 realizations. Error bars are very 
smaU. Food webs (x), C. Elegans (•), WTW (▼), cortical 
networks (a). 



relation conditions considered in this paper. The equa- 
tions are valid even for small networks such as the cat 
cortical network {N = 53, density of links d « 0.3) de- 
spite the fact that formulas are derived in the thcrmo- 
dynamical limit. Error bars in most cases are very small 
and only food webs exhibit some larger fluctuations. 



V. SUMMARY AND DISCUSSION 

In summary, we have studied the influence of 1-node 
and 2-node degree correlations on the reciprocity of net- 
works with arbitrary degree sequence. We find that, for 
a large class of complex networks, correlations account 
for a very large part of the observed reciprocity, explain- 
ing it almost completely in typical cases. In general, the 
contribution of correlations to r is nontrivial and largely 
depends on the type of correlations involved, revealing 
a hierarchy of correlation classes that contribute to r in 
different levels of approximation. As observed, this in- 
fluence can span over orders of magnitude in some real 
networks. Our analytical estimations arc proved as valid 
expectation values of reciprocity by comparison to en- 
semble averages of random networks which preserve de- 
sired types of correlations. Both from a theoretical and 
a practical point of view, it would be highly desirable 
to extend the current work and obtain expected values 



of other network measures, e.g. clustering cocfflcient, 
average pathlength, motif profiles, etc. following a sim- 
ilar philosophy: analytical expressions should be com- 
putable using only information that can be directly mea- 
sured from the specific real network under study. In the 
absence of such results, for the moment, the numerical 
methods introduced in this paper arc a useful tools for 
significance testing of network measures under conditions 
of prescribed degree sequences and degree correlations. 

The dynamical influence of degree correlations has 
been studied for epidemic spreading [23| and for synchro- 
nization in undirected networks [24| . In the light of our 
results, which clearly relate reciprocity and degree corre- 
lations, it is important to also investigate the influence 
of reciprocity on those phenomena for the broader class 
of directed networks. 

In order to perform satisfactory modeling, the key pa- 
rameters governing network growth and evolution need 
to be identified. Therefore, theoretical understanding of 
the interplay between different topological properties is 
necessary to distinguish between the significant measured 
values and those expected as by-products of other prop- 
erties. In this paper, we have explored the interplay be- 
tween degree correlations and network reciprocity. Pre- 
vious efforts in this direction include relations between 
degree correlations and clustering coefficient [11] . Degree 
correlations are also expected to be relevant for other net- 
work measures, e.g., pathlength, network motifs, modu- 
larity, etc. Therefore, further work is desirable to extend 
the analytical and experimental approaches presented in 
this paper. No doubt, quantification of similar struc- 
tural interrelations will significantly elucidate the essence 
of the structure-function-evolution interplay in complex 
networks. 
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